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We compute the (zero frequency) Hall coefficient Rh, and 
the high frequency Hall constant R* H for the strong coupling 
Hubbard model away from half- filling, in the d = oo/ local ap- 
proximation, using the new iterated perturbation scheme pro- 
posed by Kajueter and Kotliar. We find quantitative agree- 
ment of our R* H with the QMC results obtained in two di- 
mensions by Assaad and Imada [Phys. Rev. Lett. 74, 3868 
(1995)]. However while the sign of Rh is quite accurately re- 
produced by R* H the doping dependence of its magnitude at 
large U is not. We report results for the complete tempera- 
ture (T), doping (x) and U dependence of Rh and R* H and 
discuss their possible relevance to doped cuprates. 

PACS numbers: 71.28.+d,71.30.+h,72.10-d 



Hall measurements are normally regarded as a stan- 
dard probe for the density and sign of charge carriers 
in electronic systems. Within a Boltzmann, continuum 
(eg ~ k 2 ) picture, the Hall coefficient, Rh, is given by 
1/ne, where n is the density and e the charge of the 
carriers. In the presence of a lattice, and for weak scat- 
tering, Rh depends on the curvature (d 2 ej:/dk 2 , etc) of 
the Fermi surface (FS) jij and need not have the same 
sign as e; an " electron- like" (positive curvature) FS im- 
plies Rh < and a "hole-like" (negative curvature) FS 
implies Rh > 0. 

An obvious question of great interest is whether this 
correspondence would survive for strong correlations. 
Experiments on some of the high T c systems [^],[| reveal 
a hole-like FS Q while Rh changes from electron-like to 
hole-like with electron doping ||. In addition, both the 
doping and temperature dependence of Rh are "anoma- 
lous" : Rh shows unbounded growth as the " half-filled" 
state is approached, while the temperature dependence 
is roughly 1/T, with a saturation as T — ► T c , the su- 
perconducting transition temperature, from above. The 
physics of the cuprates is believed to be governed by 
strong electron-electron repulsion (large U, in the Hub- 
bard context) and low dimensionality [d = 2). This has 
motivated the search for controlled approximations that 
help to isolate the generic correlation affected features in 
transport in strongly correlated models. For the Hubbard 
model in d = 2 the only controlled calculations that exist 
are for the high frequency Hall constant, R* H , defined by 



R H = I™ lim v x , v (uj)/[BctI (uj) 

W— >oo 



(i) 



Unlike the usually measured (zero frequency) Hall co- 
efficient, Rh, which involves low energy scattering pro- 
cesses, R* H defined above involves only the reactive re- 
sponse. Nevertheless, as demonstrated first by Shastry 
et al || , R* H encodes information about correlation ef- 
fects in the system and (unlike Rh) is accessible within 
controlled calculations. For the Hubbard model it has 
been computed through Quantum Monte Carlo and 
finite cluster ^ calculations, and for the t — J model 
through high temperature series expansion ||, and clus- 
ter calculations ||. Finite size effects and analytic con- 
tinuation problems which hinder low frequency transport 
calculations in these contexts do not significantly affect 
the results on R* H . 

One significant result of these studies is the finite tem- 
perature crossover from electron-like to hole-like R* H as 
the hole doping concentration, x, away from half-filling 
is reduced below ~ 30%. This is consistent with some 
of the high T c data on Rh an d raises the intriguing 
possiblity that R* H may mimic all the features of Rh- 
Hence it would be interesting to calculate both Rh and 
R* H for a strongly correlated model within a single scheme 
and clarify the connection between them. The recently 
developed d = oo or local approximation (^Jl^l does pro- 
vide us with a controlled scheme for studying this con- 
nection in a strongly interacting system. We report and 
discuss the results of such a study of the Rh and R* H for 
the doped Mott insulator phase of the Hubbard model 
in this paper. Our primary results can be summarised 
as follows, (i) The results for R* H computed in the local 
approximation are quantitatively consistent with QMC 
results j?J obtained for d=2. (ii) There is a complete cor- 
respondence of the sign of Rh and R* H computed within 
this approximation at intermediate T over a wide range 
of x and U. (Hi) However, Rh and R* H , differ in sign as 
T — > 0, and their magnitudes have qualitatively different 
dependences on x, for intermediate T and large U , close 
to half-filling. 

Our calculations are based on the (by now well known) 
mapping of Hubbard model in the d — > oo limit to a self- 
consistently embedded impurity problem PJlCj], At half- 
filling this model shows a transition from a Fermi- 
liquid to a Mott insulating phase at U — U c ~ 2.5 (in 
units of D, the half-bandwidth, which we set to 1) at 
T = 0, as long as there is sufficient frustration in the lat- 
tice to suppress magnetic order. Away from half-filling 
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the system is metallic, and there is no evidence of mag- 
netic order Jl2]| , Most of this has been learnt Jl3| using 
a combination of Quantum Monte Carlo, exact diagonal- 
isation, the NCA [O], and iterated perturbation theory 
(IPT) |I(^[l5],[l6| schemes for solving the impurity prob- 
lem. Of these, the IPT scheme, because of its simplic- 
ity, has been particularly useful in clarifying the physics 
at half-filling. It reproduces Fermi-liquid behaviour for 
U < U c and is exact in both the atomic (t — 0) and 
the band (U = 0) limits. Although it does not repro- 
duce the exact critical behaviour as U — > U c at T = 0, 
it gives excellent results at intermediate temperatures 
and allows one to directly calculate the continuous real- 
frequency spectra of the system. The absence of a similar 
approximation away from half-filling, where the proxim- 
ity to a correlation driven insulating state is expected 
to profoundly affect the "metallic" behaviour, had until 
recently been a bottleneck. A recent paper by Kajueter 
and Kotliar JEtJ achieves a significant advance by devis- 
ing an interpolation scheme which is exact in the atomic 
and band limits at any filling and obtains the correct low 
energy physics by imposing the Friedel sum rule. Our 
calculations are based on this new scheme, which as orig- 
inally proposed was for T = 0, but has been adapted by 
us to T in a straightforward way. 



We refer the reader to the original paper |17| for the 
detailed formulation of this scheme and comparison with 
exact results. For our purposes it is enough to note that 
in this scheme the self-energy at T = is approximated 
as T, int (w) — AS 2 ( W )/[1 — BT, 2 (uj)] where S 2 is the ana- 
lytically continued second order self-energy 
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and 



have been chosen so as to reproduce 



the atomic limit and the high frequency behaviour of 
the self-energy at any filling; n = /_„ pc(^)duj and 

n Q = §°_ oo pg(uj)duj, where Pg( w ) an d Pg{^>) are the 
spectral functions of the full local propagator G(uj) and 
the "bare" local propagator G(<jj) respectively; = 
w + A*o — A(oj), where A(w) is the hybridisation of the 
self-consistently embedded impurity problem and po is a 
free parameter fixed by imposing the Friedel sum rule. 
At T = and for a given U, the chemical potential p, 
" quasi chemical -potential" po and occupation number n 
are related uniquely through the interpolation and sum 
rule constraints. 

The ground state averages n and no which enter the in- 
terpolating self-energy T, int (ui) trivially generalise to T ^ 
as thermal averages (n = (h) = J^° oo pQ(uj)f(ui)dLJ, 
f = [exp{ijj/T) + l] -1 ;, etc ) and still yield a sensible 
interpolation. For T ^ and a desired n we fix po(n,T) 



as po(n,T = 0) and solve for p(n,T) self-consistently. 
Thus once the connection between p, n and po is estab- 
lished at T = for a given U, the T / calculations 
only involve the d = oo consistency loop and are compu- 
tationally quite affordable. We compute the self-energy, 
of the lattice problem as the self-consistent self- 
energy Ej nt (u;) of the impurity problem directly for real 
frequencies, and use it to study Rh and R* H for large U, 
(U/D ~ 4 — 6) as a function of T and x. Rh is com- 
puted from the conductivity tensor a a p{uj = 0). In the 
d — > oo /local approximation the transport coefficients do 
not involve vertex corrections and the dc conductiv- 
ity is given by 

a xx {0) = c xx fdep (e) f dioA 2 {e,Lu)(~df /da) (3) 



where c xx = e 2 "ir / (dhao) ,ao is the lattice spacing, po(e) = 
(2/nD)y/l - {e/D) 2 is the bare density of states (DOS), 
and A(e, u>) = —-^Im{u + p — e — E(w)] -1 is the spectral 
function. The transverse conductivity is given by 



a xy (0) = Bc xy J dep {e)e J dcuA 3 {e,uj)(-df /du) (4) 



where c x 



e\ -K 2 ao/ C&d 2 Ti z ). Rh is then given by 
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a xy (0) /[Ba 2 x (0)]. For R* H we use the expression for a 
hypercubic lattice M 



R*H = ( a o/M)y^ cosk x cosk y (nk,a)/\y cosk x 



\ n k.a / 



(5) 



In the spirit of the local approximation we compute 
(n/c.o-), the momentum distribution function, from our 
self-energy and do the k sums over a 2d Brillouin zone. 

We first discuss R* H (x, T) since 2d QMC and cluster 
calculations || provide us with results to compare with. 
For large U, (U/t ~ 16) QMC results are available only 
above T/t ~ 0.25, (T - 0.06 in our units, D - 4t) and 
there is a clear peak structure in R* H (T) with the location 
of the maxima T^ x c gradually shifting to lower T as U 
increases at a fixed doping (x — 0.05). This is consistent 
with our results; but as T — > we find an "anomalous" 
region (R* H > 0) in the U — x plane which could not 
be accessed within QMC but is clearly seen in cluster 
calculations ||. The detailed comparison is shown in 
Fig.l. 

Exploring the detailed x, T dependence of R* H and Rh 
for U/D = 4 (roughly U = I6t in 2d), we find that R* H re- 
mains hole-like down to T = for x < 0.12 (Fig. 2a) while 
for x > 0.20, R* H is electronlike for all T. Rh however is 
electron-like for all x at T = (Fig. 2b). This, as we dis- 
cuss below, is a consequence of the Friedel sum rule which 
forces Rh to preserve its non-interacting, band structure 
determined, value at T = 0. For small x however there 
is a rapid crossover to a hole-like Rh with increasing 
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T. Fig. 3 shows the crossover temperature T cr as located 
from the numerically evaluated Rh ■ T cr ~ x 2 for small, x 
(Fig. 3), and rises more steeply clS X ^ Xr rsJ 0.25, beyond 
which there is no crossover. For x < x c ~ 0.25 there is a 
high temperature branch of T cr (x) across which there is 
a "re-entrant" change of sign in Rh- For x > 0.25 Rh, 
like i?^, is electron-like for all T. While the T depen- 
dence illustrates the similarity in the sign of Rh and R* H 
at intermediate T (T > 0.15 ~ O(J) in 2d) the T -> re- 
sult illustrates their essential difference. This difference 
arises from the fact that, for T — ► 0, is a Fermi sur- 
face object, completely dominated by the coherent part of 
the spectra (which leads to the divergent <J a p), while R* H 
being a "high-frequency" object gets contribution from 
all over the Brillouin zone. 

The fact that at T = all correlation effects cancel out 
of Rh @ is easily seen by using the uj — > 0, T — > 
Fermi liquid parametrisation of the spectral function : 
A(e,uj) ~ T/((m*uj + p, - e) 2 + 7r 2 r 2 ) where T = 
-(1/tt)Ito£(cj = 0), m* = 1 - (dE R /du) u = and 
/2 = - S fi (.x,T,w = 0). Since -df/duj -► 5{lu) 

as T -> 0, er^ ~ p (fi)/T, a xy ~ Bjlp (fi)/T 2 and 
i?# ~ /l so long as r << 13. Since the sum rule requires 
/2(x, T = 0) = /i(x, T = 0,U = 0) = fJ,o(x), Rh is pinned 
to its band value at T = 0. The small x behaviour of the 
lower T cr (Fig. 3) can also be understood within this phe- 
nomenology. For x — > (accessible within our numerics 
so far upto x ~ 0.06) p(x,T) increases rapidly with T. 
We find p(x,T) ~ Ho{%) + aT 2 /x 3 , where a > depends 
on [/. Since the free band chemical potential po(x) ~ —a; 
for x — > 0, the crossover scale T cr ~ a; 2 . The rapid rise of 
/2 with increasing T is primarily due to the reduction of 
Efl(cc, T,u = 0). This in turn is due to the redistribution 
of spectral weight in ImS(w) which rises rapidly at low 
frequencies as T increases, causing the disappearance of 
the resonant feature in pc{u) at the Fermi level. 

The difference between Rh and R* H is further high- 
lighted in their doping dependence as x — > at inter- 
mediate T. Rh ~ g(T)/x with g(T) increasing with T. 
This behaviour is seen down to the lowest accessible x 
(~ 0(0.04) at these T), despite particle-hole symmetry 
which should force Rh{T) = at x — 0. This indicates 
that due to the large finite T resistivity in the "metal" 
for x — > (and the absence of any coherence/sum rules 
unlike at T = Q), Rh "diverges", as if the doped holes 
could be treated as effective carriers within a " Boltzmann 
like" scenario. R* H on the other hand shows a peak as x 
reduces and then a clear turnaround, with R* H ~ x for 
x — > 0. Thus as the "Mott" phase is approached for 
T 7^ 0, i?^ is "divergent" (with possibly an exponen- 
tially small turnaround scale because a a p is finite) while 
R* H is benign. This feature of R* H has been seen within 
QMC § and confir ms that for a one band model with 
particle- hole symmetry R* H is regular as x — > 0. The "di- 
vergence" in Rh as x — > had been seen earlier within 



slave-boson gauge theories based on the f — J model [^3| , 
though they yield only a hole- like Rh (i-e no sign change 
even as x — * 1, i.e n — ► 0). Within the d — ► oo scheme 
Pruschke e£. a/ pl| have computed Rh using a combi- 
nation of QMC and the NCA. Our results are consistent 
with the QMC data while avoiding the problems of an- 
alytic continuation and access low T, large U regimes 
much better. 

It is natural to compare our results with the Rh data 
on the cuprates (2|||| since the cuprates are believed to 
be doped "Mott" insulators. Although certain features 
of the data, such as the sign change with increasing x || , 
" re-entrance" 0| , and the 1 jx behaviour of Rh |||| are 
recovered within the local approximation the T depen- 
dence of Rh, and the characteristic temperature scales 
are not captured. Even taking into consideration that 
the T = metallic state is not accessible due to super- 
conductivity, experimentally Rh decreases (roughly as 
l/T) beyond a T scale which is O(50K - WOK). In 
our calculation Rh decreases with increasing T only be- 
yond a temperature scale ~ 0(D /10). As demonstrated 
by its success in understanding metal-insulator transi- 
tion in 3d transition oxides p5| the local approximation 
quantitatively captures the physics of strong local cor- 
relations. To the extent that one believes that the low 
energy physics of the cuprates can be described by the 
large U Hubbard model, one must attribute deviation 
of the measured from the computed Rh to one or more 
of the following features which are left out in our treat- 
ment: (a) bandstructure effects, including particle-hole 
asymmetry and the proximity to the Fermi level of the 
van Hove singularity in the quasi 2d DOS. (b) the k de- 
pendence of single particle self-energy, and (c) vertex cor- 
rections to the conductivity, which do not vanish when 
(b) is relevant. What is nice about the local approxi- 
mation is that it has allowed a clean separation of the 
local correlation effects on Rh, which we have discussed 
in this paper, from the effects attributable to these other 
complications, which we hope to study in future work. 

We thank P. N. Mahesh, B. S. Shastry and T. V. Ra- 
makrishnan for discussions, and Alexander Punnoose and 
Prabhakar Pradhan for help with the figures. 
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FIG. 1 R* H in units of ag/|e|, for different U/D. Trian- 
gle U/D = 4; inverted triangle U/D = 3; cross U/D = 2. 
Inset: QMC data 0, same legend for U/At, 

FIG. 2 R* H (a) and R H (b) in units of ajj/|e|, for U = 4. 

FIG. 3 Crossover temperature T cr (x) for U = 4. The 
line is a quadratic fit for x — > 0. "normal" and "anoma- 
lous" refer to Rh < and Rh > respectively. 

FIG. 4 R H (a) and R* H (6) in units of ag/|e| for U = 4. 
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